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Abstract 

We study the localization properties of the eigenstates in the Kronig- 
Penney model with weak compositional and structural disorder. The 
main result is an expression for the localization length that is valid 
for any kind of self- and inter-correlations of the two types of disor- 
der. We show that the interplay between compositional and structural 
disorder can result in anomalous localization. 

Pacs numbers: 73.20.Jc, 73.20.Fz, 71.23.An 

Recently, much attention was paid to low- dimensional disordered mod- 
els with long-range correlations in random potentials. Apart from the the- 
oretical aspects, the interest on this issue has increased significantly due 
to the possibility of constructing random potentials with specific correla- 
tions which result in a strong enhancement or reduction of the localization 
length [H [21 El m El [6] . These new effects allow for the fabrication of electron 
and optic/electromagnetic devices with desired anomalous transport proper- 
ties. As was shown analytically [21 El HJ and confirmed experimentally [5l E] , 



1 



one can arrange prescribed windows of energy with perfect transmission (or 
reflection) of scattering waves. 

One of the most important models, both from the theoretical and exper- 
imental point of view, is the Kronig-Penney (KP) model, which was intro- 
duced long ago to analyze electronic states in crystals [7]. Since the '80s, this 
model has attracted considerable attention because it provides a convenient 
description of superlattices (see, e.g., [S] and references therein). Modifica- 
tions of the standard Kronig-Penney model have been suggested for a study 
of the physics of random and quasi-periodic systems with various applica- 
tions, see, e.g., [9]. Recently, the Kronig-Penney model has been used to 
discuss the possibility of selective transmission in waveguides (see E] and 
references therein). 

In this paper we study the KP model with two types of weak disorder. 
Disorder of the first kind, or "compositional", is due to small variations in 
strength of the delta-shaped barriers. In addition, the spacings between the 
barriers can be also randomly perturbed (the so-called "structural" disorder). 
Our interest lies in the interplay of these two kinds of disorder which can 
exhibit both self-correlations and mutal correlations. Our goal is to derive a 
formula for the localization length, and to analyze it. 

The stationary Schrodinger equation for the eigenstates ip{x) has the form 

-— ^"(X)+ J2 Un6{x-Xn)Hx)=E^lj{x), (l) 
n=-oo 

where Un = U + Un and amplitude and position of the n-th. 5-barrier. 

In what follows we use units in which /2m = 1; we can thus write the 
energy of the eigenstates a.s E = where q is the electron wavenumber. 

The positions of the 5-barriers are assumed to be slightly shifted with 
respect to the lattice sites, Xn = na + an, where a is the lattice step. The 
variables Un represent fluctuations of the barrier strength around the mean 
value U . Our analysis is restricted to the case of weak disorder for which 
both variables Un and a„ have zero average, = and (a„) = 0, and 
small variances, q'^{aD -C 1 and {uD -C U^. We remark that the condition 
q'^{aD ^ 1 implies that the energy must be low on a scale set by l/(a^). 
Note that in contrast with many previous studies, both variables are random 
and may have stationary correlations: our main interest lies in how these 
correlations shape the properties of the localization length /loc(-E) of the 
eigenstates. 
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It is convenient to introduce the relative displacements of the barriers, 
= 0'n+i~0'n^ having zero mean, (A„) = 0, and small variance, q'^(A^) <^ 1. 

Apart from the first two moments of the random variables and m„, one 

has to give the binary correlators, 

Xi{k) = {unUn+k)l{ul) 

X2{k) = (A„A„+,)/(A2) (2) 
Xsik) = (u„A„+fc)/(M„A„). 

We will not attribute specific forms to the correlators Xi{k); we simply as- 
sume that they depend only on the index difference k because of the spatial 
homogeneity in the mean of the model and that they are even functions of 
k. 

It is worthwhile to note that Eq. ([T]) can be treated as the wave equa- 
tion for electromagnetic waves in a one-dimensional (ID) waveguide with 
wavenumber q = u/c. Therefore our results are equally applicable to the 
classical scattering in optical and microwave devices of the Kronig-Penney 
type with correlated disorder. Our model is also equivalent to a classical os- 
cillator with a parametric perturbation constituted by a succession of 5-kicks 
whose amplitudes and time-dependence are determined by f/„ and a„. This 
correspondence allows one to cast Eq. ([T]) in the form 

oo 

x+[q^- E UJ{t-Q]x = 0. (3) 

n=— CO 

Our analysis is based on the Hamiltonian approach [TOl [11] according to 
which the spatial structure of eigenstates of the KP-model can be analyzed 
by exploring the time evolution of the kicked oscillator described by the dy- 
namical equation Such a dynamical approach considers the Schrodinger 
equation as an initial- value problem and can be treated as a modification of 
the transfer matrix approach. 

Integrating the dynamical equation ([3]) between two successive kicks, one 
obtains the map 

Xn+i = [{Un/ q) sin (/i + + cos (/i + /i„)] x„ 

+ (l/g)sin(/i + /i„)p„ 
Pn+i = [Un COS (/i + /i„) - g sin {jj, + Xn 

+ COS (/i + Pn 

where fi = qa and fin = Q'A„, and the values Xn and p„ refer to the instant 
before the n— th kick. 
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The evolution of the dynamical map (jl]) can be analyzed as follows. First, 
we make a weak-disorder expansion of Eq. (jlj) , keeping only first- and second- 
order terms. The expansion is straightforward; the resulting equations, how- 
ever, are lengthy and we omit them here. As a second step, we perform 
a canonical transformation — > {Xn,Pn), such that the unperturbed 
motion reduces to a simple rotation in the phase space of the new vari- 
ables [3j. Such a trick allows one to eliminate the effect of the periodic kicks 
with constant amplitudes U. This can be done with the use of the canonical 
t r ansf or mat ion , 

Xn = acos(/i/2)X„ + (ga)"^sin(/i/2)P„ 
Pn = -gasin(/i/2)X„ + a~^cos(;u/2)P„ 

where the parameter a is defined by the relation 

1 sin/x- g(cos/i- 1) 

ct = . 

sin /i — 1^ (cos /i + 1) 

Note that, due to the transformation ([5]), the new variables X„ and P„ have 
the same dimension. 

In the absence of disorder, i.e., for m„ = and A„ = 0, the rotation angle 
7 between successive kicks is determined by the relation, 

cos 7 = cos u H sin u with 7 = ka. (6) 

2q 

In terms of the Kronig-Penney model, k is the Bloch wavevector, 7 is the 
phase shift of the wavefunction within the lattice step a, and Eq. (jH]) defines 
the band structure of the energy spectrum. 

It should be pointed out that the transformation ([5]) is well-defined for 
all values of the rotation angle 7 other than 7 = and 7 = ±7r for which a 
either vanishes or diverges. In other words, our approach fails at the center 
and at the edges of the first Brillouin zone, i.e., at the edges of the allowed 
energy bands of the KP model. However, the approach works well in every 
neighborhood of these critical points. 

To proceed further, it is useful to pass to the action-angle variables 
{Jni(^n)i with the transformation 

X„ = y2Xsin6'„, Pn = yfn^cosOn 
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and to represent the Hamiltonian map (jlj) in terms of the new variables. 
Leaving aside mathematical details, we give here the final expression, 



Jn+l 



J 



where 



Dl 



= ^„ + 7 - Hi - cos (2^„ + 7)] Un 
+ i - COS (2e„ + 27)] A„ 

1 + sin (26'„ + 7) u„ - sin (26'„ + 27) A, 



(7) 



+ |[l-t;cos(2e„ + 27)] A2 



+ 



(8) 



cos {26n + 7)] u. 
- [cos 7 - COS {29n + 27)] M„ A„. 

Here, v = [ga^ + (ga^)""*^] gsin7/f/, and we have introduced the rescaled 
random variables, 



Ur. 



sm /i 
gsin7 



and Ar 



U 



sm7 



■A, 



In Eqs. ([7]) and ([8]) we have kept only the terms of the weak-disorder ex- 
pansion which are necessary to compute the localization length within the 
second-order approximation. We remark that the angle variable evolves in- 
dependently of the action variable. 

The inverse localization length for the KP model can be computed as 
the Lyapunov exponent A, 



'loc 



1 ^ 

A = lim In 

^-00 Na 





= (-\n 











which, in terms of the dynamical map ([7]), can be written as [TT 



A 



\2a 



J n. 



(9) 



By expanding the logarithm of D"^, one gets 

A = ^({sin(2e„ + 7)u„-sin(2^„ + 27) A„ 

+ l[l-2v cos (20„ + 27) + cos (AOn + 47)] Al 
+ i[l-2cos(2^„ + 7) + cos(4^„ + 27)]M2 

2 cos (2ft„ + 27) + cos {AOn + 37)] UnAn} ). 



(10) 
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Now, in order to obtain the Lyapunov exponent A, we have to perform 
the average over the phase 6n and the random variables Un and A„. To the 
second order of perturbation theory, one can neglect the correlations between 
9n and the quadratic terms u^, A^, and u„A„. Hence for the summands in 
Eq. ffTOj) which contain these quadratic terms, one can compute separately 
the averages over 6'„ and over the random variables ti„ and A„. 

In analogy with the Anderson model (see details and references in 
it can be shown that for our purposes it is safe to assume that the invariant 
measure of the phase is a flat distribution, p{9) = l/{2n). The assumption 
holds for all values of 7, except for 7 = ±7r/2 where a small modulation of 
the invariant measure results in an anomaly for the localization length. The 
situation with these values of 7 is similar to that known for the standard ID 
Anderson model at the center of the energy band, and the correct expression 
for A can be obtained following the approach of [11] . 

It should be stressed that weak modulations of p{9) arise also for other 
"resonant" values, 7 = mn/r, with m and r integers prime with each other 
and r > 2. However, these modulations do not influence the value of A, 
because the expression to be averaged in Eq. flTU]) has no harmonics higher 
than 46. Thus, our further analysis is valid for all values of 7 except 7 = 
and 7 = ±71 (i.e., the edges of the energy bands) and 7 = ±tt/2. 

After averaging, the expression for the Lyapunov exponent takes the form. 



A 



, (^n) + (An)-2(n„A„)cOS7 
f 1 ~ 

+ — (M„sin(2^„ + 7)) - — (A„sin(2ft„ + 27)). 



(11) 



In order to compute the noise-angle correlators in Eq. f[TT|) . we generalize 
the method used in [2]. Specifically, we introduce the correlators = 
(ttn exp (i26'„_fc)) and Sk = {Anexp {i29n-k)) ■ Both correlators satisfy re- 
cursive relations that can be obtained by substituting the angular map of 
Eq. ([7]) into the definitions of Vk-i and Sk-i- The recursive relations allow 
one to obtain the correlators tq and Sq, whose imaginary parts represent 
the noise- angle correlators in Eq. flTTl) . As a result, we arrive at the final 
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expression for the Lyapunov exponent, 



A 




{ul)Wi + {Kl)W2 - 2(m„A„) cos7l^3 



(12) 



where the functions 



oo 



Wi = Wi (27) = 1 + 2 ^ Xi{k) cos(27fc) {i = 1, 2, 3) 



(13) 



k=l 



are the 27— harmonics of the Fourier transform of the binary correlators 
Xi{k)^ see Eq. ([2]). We should stress that Eq. f[T^ has been derived without 
invoking the Born approximation, E it describes the tunneling regime 

for E < U as well as the scattering one, E > U. The only constraint is the 
weakness of both types of disorder, {uD <^ f/^ and g^(A^) <^ a^. 

Let us first discuss the structure of expression fll2p for the case in which 
there are no correlations between Un and A„. It is quite instructive that 
for weak scattering, U <^ q, (therefore, 'j ^ fi) with purely compositional 
disorder, (i.e., A„ = 0), the Lyapunov exponent takes the form A ~ g^^i) 
(see j6jj). This is equivalent to the well-known result for weak scattering in 
continuous ID potentials, A = ^'^W{2q), where {V'^) is the variance of the 
random potential and W{2q) is the 2g— component of the power spectrum of 
the potential. 

In the other limit case of structural disorder (i.e., m„ = 0), the expression 
for the Lyapunov exponent, A = gQ^n^^ (^n)^2; was obtained in [3]. One 
can see that it is similar to the expression for the tight-binding Anderson 
model with diagonal disorder, A = ^J^^.J^ ^ W{2n), where (e^) is the variance 
of the random site-potential, fi is the phase shift of the wavefunction between 
two sites -related to the energy by the dispersion relation E = 2 cos /i-, and 
W{2^) has the same meaning as Wi. 

Expression f|T2l) for the inverse localization length allows one to estimate 
the relative importance of the structural and compositional disorder for the 
transport properties of a finite sample. As can be seen, spatial correlations 
of the variables Un and A„ can enhance or suppress the localization length 
in comparison with the case of uncorrelated disorder. 
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Specific long-range correlations can make the Fourier transforms fll3l) van- 
ish in prescribed energy windows, so that the Lyapunov exponent (IT^ also 
vanishes in the same energy intervals. A method for the construction of ran- 
dom potentials with given binary correlators Xi(^) was described and tested 
in [21 Hj. Following the same approach, one can use formula f|T2|) clS db stjcirtj" 
ing point for the fabrication of devices with prescribed anomalous transport 
characteristics. 

It is interesting to relate the properties of the KP model with known 
results for ID tight-binding models with both diagonal and off-diagonal dis- 
order. It can be shown that, after eliminating the momenta pn from the 
map (jl]), one obtains the equation 



[1/ sin(/i + fin)] ipn+i + [1/ sin(/i + tp. 



n-l 



COt{n + /i„) + COt{n + /i„_i) + i (f/ + Un) 



(14) 



where V'n = Xn- For weak disorder one can expand the coefficients of Eq. ( 031) 
in powers of gA„ and get 



1 + 1/2 + cot^ /il (A^ ) - gA„ cot /i] V^„+i 
1 + (1/2 + cot^ /i) (A^ ) - gA„_i cot /il ^l). 



•->n-\ 



2cos/i + f/^ + 2g2£^^|ii(A2) 

-Un '^n- 



sm ^ 



(A„ + A„„i) + 



(15) 



One can see that the right-hand side of this relation vanishes under the 
conditions 



U 



Ur 



2q cos fi 



sm fi 



1 + 



sin^ /i 



sin^ /i 



(A„ + A, 



n-l) 



(16) 
(17) 



The first condition determines the energy as a function of the mean field 
U and of the variance (A^). The second condition establishes non-trivial 
correlations between the compositional and structural disorders. Under these 
conditions Eq. (fT5|) reduces to 



1 + q^ (l/2 + cot' fi) (Al) - gA„ cot fi^ , 
^1 + q^ (l/2 + cot' /i) ( A' ) - gA„_i cot /i] ^„_i = 



(18) 
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and this identity can be written as the Schrodinger equation 

(7 + 7„+i) tpn+l + (7 + In) = (19) 

for the Anderson model with purely off-diagonal disorder at the center of the 
energy band. 

As is known , the model fll9p exhibits anomalous localization. Specif- 
ically, the band-center electronic state is localized but decays away from the 
localization center no as ipn ~ exp (^—A^J\n — no| j where A is some constant. 
Thus, the interplay between compositional and structural disorder can give 
rise to anomalous localization in the KP-model. 

We have to note that the above conclusion about the anomalous localiza- 
tion cannot be drawn directly from the general formula fll2p for two reasons. 
First, the zero- value of the Lyapunov exponent leaves open the question of 
whether the corresponding electronic state is extended or anomalously lo- 
calized. Second, the conditions f|T6l) and f|T71) imply that the Bloch vector 
k lies in a neighborhood of the points ±7r/2a, where our derivation of the 
Lyapunov exponent may be invalid. Indeed, when conditions ( fT6l) and ( ITTI) 
are met, the requirement of weak compositional disorder, (m^) -C U"^ , leads 
to g^(A^) <C sin^ /i . Taking into account this inequality and the relation ([6]) 
with U given by (fT6!) . one can see that the Bloch wavenumber is actually 
close to the resonant values ±7r/2a. 

In conclusion, we have derived the expression for the localization length 
in the KP model for the general case when both the amplitudes and the 
spacings of the 5-barriers are random variables. Our consideration takes into 
account correlations of each type of disorder with itself, as well as correlations 
between the two disorders. The obtained expression depends on the binary 
correlators only and thus opens the way to the construction of random poten- 
tials with desired characteristics of the electron (or electromagnetic waves) 
transmission through finite samples. The most important application of ex- 
pression (fT2!) lies in its use for the fabrication of devices with prescribed 
energy windows with perfect transmission (or reflection) produced by long- 
range correlations of the disorder, an effect recently observed experimentally 
in microwave waveguides [51 We have also found that for specific cor- 
relations between the two kinds of disorder, the Kronig-Penney model has 
anomalously localized eigenstates, similarly to the Anderson model with off- 
diagonal disorder. 

The authors gratefully acknowledge the financial support of the Coecyt 
grant CB0702201-X. 
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